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Bilayer quantum Hall system (BLQH) differ from its single layer counterparts (SLQH) by its 
symmetry breaking ground state and associated neutral gapless mode in the pseudo-spin sector. Due 
to the gapless mode, qualitatively good groundstate and low energy excited state wavefunctions at 
any finite distance is still unknown. We investigate this important open problem by the Composite 
Boson (CB) theory developed by one of the authors to study BLQH systematically. We derive the 
ground state, quasi-hole and a pair of quasihole wavefunctions from the CB theory and its dual 
action. We find that the ground state wavefunction is the product of two parts, one in the charge 
sector which is the well known Halperin's (111) wavefunction and the other in the spin sector which is 
non-trivial at any finite d due to the gapless mode. So the total groundstate wavefunction differs from 
the well known (111) wavefunction at any finite d. In addition to commonly known multiplicative 
factors, the quasi-hole and a pair of quasi-holes wavefunctions also contain non-trivial normalization 
factors multiplying the correct ground state wavefunction. We expect that the quasi-hole and 
pair wave function not only has logarithmically divergent energy, well localized charge distribution, 
but also correct interlayer correlations. All the distance dependencies in all the wavefunctions are 
encoded in the spin part of the ground state wavefunction. The instability encoded in the spin part of 
the groundstate wavefunction leads to the pseudo-spin density wave formation proposed previously 
by one of the authors. Some subtleties related to the Lowest Landau Level (LLL) projection of the 
CB theory are also noted. 

I. INTRODUCTION 

The Wavefunction approach has been very successfully applied to study single layer quantum Hall ( SLQH ) 
systems at Laughlin series v — 2f } +1 [1] and Jain's series at v = 2s ^ ±1 [2]. One of main reasons for the success of the 
wavefunction approach in SLQH is that there is a gap in the bulk, suitable wavefunctions [1,2] can describe both the 
groundstate and low energy excitations quite accurately. Its accuracy can be checked easily by exact diagonalization 
in a finite size system whose size need only go beyond a few magnetic length. Spherical geometry can be used to get rid 
of edge state effects quite efficiently. In general, trial wave function approach is a very robust approach to study SLQH 
and multi-components systems as long as there is a gap in the bulk. The gap protects the many properties of the 
system such as charge density distributions and energies from being sensitive to some subtle details of wavefunctions. 

However, the situation could be completely different in Spin-polarized Bilayer Quantum Hall system ( BLQH ) at 
total filling factor vt = 1. This system has been under enormous experimental and theoretical investigations over the 
last decade [3]. When the interlayer separation d is sufficiently large, the bilayer system decouples into two separate 
compressible v — 1/2 layers [4]. However, when d is smaller than a critical distance d c \, even in the absence of 
interlayer tunneling, the system undergoes a quantum phase transition into a novel spontaneous interlayer coherent 
incompressible phase which is an excitonic superfiuid state ( ESF ) in the pseudospin channel [5-7]. In [10], Halperin 
proposed the (111) wavefunction to describe the groundstate of the ESF state. Starting from the (111) wavefunction, 
using various methods, several authors [11] discovered a Neutral Gapless Mode (NGM) with linear dispersion relation 
u> ~ vk and that there is a finite temperature Kosterlitz-Thouless (KT) phase transition associated with this NGM. 
By treating the two layer indices as two pseudo-spin indices, Girvin, Macdonald and collaborators mapped the bilayer 
system into a Easy Plane Quantum Ferromagnet (EPQFM) [3,12] ( which is equivalent to the ESF ) and explored 
many rich and interesting physical phenomena in this system. 

As first pointed out in [12], the (111) wavefunction may not be qualitatively good at finite d, because (111) is 
a broken symmetry state in a direction in XY plane of isotropic ferromagnet at d = instead of a easy-plane 
ferromagnet at finite d. The NGM is a hallmark of the interlayer coherent Quantum Hall state. Its existence is 
expected to dramatically alter the properties of the wavefunctions of the ground state, quasi-hole and quasi-particle. 
In [13], G. S. Jeon and one of the authors studied properties of essentially all the known trial wavefunctions of ground 
state and excitations in bilayer quantum Hall systems at the total filling factor vt = 1. The results indicated that 
qualitatively good trial wave functions for the ground state and the excitations of the interlayer coherent bilayer 
quantum Hall system at finite d are still not available and searching for them remains an important open problem. 
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Specifically, they investigated the properties of the quasi-hole wave function, meron wavefunction and a pair of meron 
wavefunction built on (111) state which have superscripts "prime" in this paper: 
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where $m is the Halpcrin's (111) wavefunction: 

JVi N 2 JVi n 2 ^ 

*m(z,w) = J[(zi - z 3 )Y[( Wi - Wj) H Y[(zi - Wj )exp(-^2 ^ |z,| 2 + K| 2 ) (2) 
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where Aq = iV 2 = AT in the balanced case and z and w are the coordinates in layer 1 and layer 2 respectively. In the 
following, we suppress the exponential factor. 

These quasi-hole and meron wavefunctions differ only by "normalization factors". As shown in [13], the normaliza- 
tion factor \zi\ is accurate only at long distance limit \zi \ — > oo limit. Near the origin, the "meron" and the "quasi-hole 
' have similar behaviors. Normalization factors have been shown not to be important in single layer Quantum Hall ( 
SLQH ) systems. However, as shown in [13], they make a dramatic difference in BLQH. Although the smallest meron 
has a localized charge 1/2 and logarithmically divergent energy, the charge of the quasi-hole excitation extends over 
the whole system and its energy also diverges linearly as the area of the system size. This indicates that the quasi-hole 
wavefunction is not a good trial wavefunction for any low energy excitations. The meron wavefunction is not a good 
trial wavefunction either, because it ignores the strong interlayer correlations [13]. It was found the energy of the 
possible wavefunction of a pair of merons in Eqn.l increases quadratically ~ \z — w \ 2 instead of logarithmically 
as the separation of the pair increases. All the results achieved in [13] indicate that qualitatively good trial wave 
functions in the interlayer coherent bilayer quantum Hall system at finite d, both for ground state and excitations, arc 
still unknown and searching for them remains an important open problem, so the wavefunction approach to BLQH 
is much more difficult and far less powerful in BLQH than in SLQH. Fortunately, effective theory approaches such 
as EPQFM approach [12,3] and Composite Boson theory approach [11,12,14-16] circumvent this difficulty associated 
with the unknown wavefunction at any finite d and are very effective to bring out most of the interesting phenomena in 
the pseudo-spin sector in this system. In fact, all these effective theories start from the insights gained from Halperin's 
(111) wavefunction which is exact at d = 0. 

In a series of papers [14-16], one of the authors developed a systematic composite boson approach to study balanced 
and im-balanced Bi-Layer Quantum Hall systems in rather details. The theory puts spin and charge degree freedoms 
in the same footing, explicitly bring out the spin-charge connection and classify all the possible excitations in a 
systematic way. Then he pushed the theory further to understand novel phases and quantum phase transitions as the 
distance between the two layers is changed. He found that starting from the well studied excitonic superfluid (ESF) 
state, as distance increases, the instability driven by magneto-roton minimum collapsing at a finite wavevector in the 
pseudo-spin channel leads to the formation of a pseudo-spin density wave (PSDW) at some intermediate distances. 
He constructed a quantum Ginsburg-Landau theory to study the transition from the excitonic superfluid (ESF) to 
the PSDW and analyze in detail the properties of the PSDW. He showed that a square lattice is the favorite lattice 
and the correlated hopping of vacancies in the active and passive layers in the PSDW state leads to very large 
and temperature dependent drag observed in the experimental. In the presence of disorders, the properties of the 
PSDW are consistent with all the experimental observations [6,8] in the intermediate distances. Further experimental 
implications of the PSDW are given. Then he extended the Composite Boson theory to study slightly im-balanced 
BLQH. In the global U(l) symmetry breaking excitonic superfluid side, as the imbalance increases, the system supports 
continuously changing fractional charges. In the translational symmetry breaking PSDW side, there are two quantum 
phase transitions from the commensurate PSDW to an in-commensurate PSDW and then to the excitonic superfluid 
state. These results explained the experimental observations in [9] very nicely. The author found that the theory 
can be easily extended to study some additional interesting phenomena in trilayer quantum Hall systems [17]. It was 
concluded in [16] that field theory approaches are much more powerful in BLQH than in SLQH. 

Obviously, the CB theory circumvent this difficulty associated with the unknown wavefunction at any finite d and 
is used to achieve all these interesting and important results at two different distance regimes without knowing the 
precise wavefunctions for the ground state and excitations. It would be interesting to use the CB theory to address the 
important and outstanding problem avoided in [13] and in all the other pervious work that finding the good ground 
state and low energy excited wavefunction for BLQH at any finite d. In SLQH, the CB theory developed in [18] was 
used to re-derive the already well known Laughlin's wavefunctions for ground state and quasihole at v = ■ As 
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previously, the gap in the bulk protects the properties of the system such as charge density distributions and energies 
from being sensitive to some subtle details of wavcfunctions. Here we are facing a more difficult and interesting task: 
to derive these unknown wavefunctions at finite d. 

The rest of the paper is organized as the following. In Sec. II, in order to be self-contained, we review briefly the CB 
approach and its dual action developed in [16] which are needed to derive the wavefunctions in the following sections. 
In Sec. Ill, using the formalism presented in Sec. II, we derive the ground state wavefunction which is different from 
the (111) wavefunction at any finite d. In Sec. IV, using the dual action presented in Sec. II, we derive the quasi- hole 
wavefunction and compare it with the "quasi-hole" and "meron" wavefunction built on (111) wavefunction listed in 
Eqn.l. In Sec.V, we derive a pair of meron wavefunction with charge 1 and compare it with the " pair meron " 
wavefunction built on (111) listed in Eqn. 1. In Sec. VI, we look at the instability in the ground state wavefunction 
as distance approaches d c \. Finally, we reach conclusions in Sec. VII. Some caveats related to the Lowest Landau 
Level (LLL) projection of the wavefunctions are also pointed out. We note that there is also an alternative approach 
in [19]. 



II. COMPOSITE BOSON APPROACH AND ITS DUAL ACTION IN BLQH 

In this section, we briefly review the formalism developed in [16] which is needed to derive the wavefunctions in 
the following sections. Consider a bi-layer system with N\ ( N 2 ) electrons in left ( right ) layer and with interlayer 
distance d in the presence of magnetic field B = V x A ( Fig.l): 

H = Hq + Hi n t 

H = J d 2 xcl{xY- 2 J c a {x) 

H int = ^ J d 2 xd 2 x'5p a {x)V oll }{x - x*)8pp{x') (3) 

where electrons have bare mass m and carry charge — e, c a , a = 1, 2 are electron operators in top and bottom layers, 
Sp a (x) = <$ a (x)c a {x) — n ai a = 1,2 are normal ordered electron densities on each layer. The intralayer interactions 
are Vn = V 22 = e 2 /«", while interlayer interaction is V12 = V21 = e 2 /e^/r 2 + d 2 where e is the dielectric constant. 
Performing a singular gauge transformation [14,16]: 

<P a {x) = e *Jd 2 x'^z') P (z') Ca{ ^ (4) 

where <j>(x — x') = arg(x — x') is the angle between the vector x — x' and the horizontal axis. p(x) — c\{x)ci(x) + 
c\{x)c2{x) is the total density of the bi-layer system. Note that this transformation treats both c\ and ci on the same 
footing. This is reasonable only when the distance between the two layers is sufficiently small. It can be shown that 
<t> a {%) satisfies all the boson commutation relations. We can transform the Hamiltonian Eqn. 3 into the Lagrangian in 
Coulomb gauge: 

t, x t, A-itiJ + ^A{x)-ha{x)) 2 , 

£ = 4>1{°t - ««o)0a + 4>\{X) <f>a{x) 

+ \j d2x ' S P(x) V +(x- x')5p{x') 

+ lf d 2 x'5p-{x)V-{x - x 1 )Sp-{x') - ^-a (V x a) (5) 
2 J 2tt 

where V± = Vl1 f Vl2 and V u = V 22 = ^f,V 12 = ^! e -« d . The Chcrn-Simon gauge field is a = J dVV^x - 
x')p(x') = fd 2 ^^^p(x'). 

In Coulomb gauge, integrating out an leads to the constraint: V x a = 2-K<p\<p a . Note that if setting V- = 0, then 
the above equation is identical to a single layer with spin in the absence of Zecman term, so the Lagrangian has a 
SU(2) pseudo-spin symmetry. The V- term breaks the SU{2) symmetry into U(l) symmetry In the BLQH at finite 
d, V- > 0, so the system is in the Easy-plane limit. 

We can write the two bosons in terms of magnitude and phase 
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<Pa = VPa + S Pa e l0 « (6) 

The boson commutation relations imply that [5p a (x),9b(x)] = ihS a bS(x — x'). After absorbing the external gauge 
potential A into a, we get the Lagrangian in the Coulomb gauge: 

C = i5p+{^d T 8+ - oo) + £^l\v8+ + - *s)V0- - a] 2 + ^Sp+V+(q)Sp+ - ^-a (V x a) 

+ l -6p-d T 6- + ^.{\ve-? + l -5p-V.{q)6p- - h z Sp- (7) 

where / = Aviv 2 which is equal to 1 at the balanced case and h z — V-p~- — V-(pi — pz) plays a role like a Zeeman 
field. 

Performing the duality transformation on Eqn.7 leads to the dual action in terms of the vortex degree of freedoms 
= ^ e ^v\dyid v Q± = JJ^ 1 ± J v2 and the corresponding dual gauge fields b^: 

C d = -iwb+e^x&bt - iA+^xdJt + ntb+r+ + ^j(d a b+ - b+) 2 + i(V x b+)V+(q)(\7 x 6+) 

- iA-^ vX d v b- + mb-F- - MV x 6") + ^j{d a bo - d b~f + i(V x b~)V-^(W x b~) 
in 

- j- f {vi - ^{dfa - d Q bJ 3 ){d (j b+ - d b+) (8) 

where = A\ ± A 2 ^ are the two source fields. It is useful to stress that the dual CS term only appears in the 
charge sector. 

For simplicity, we only consider the balanced case. Putting v\ = v<i = 1/2 and h z = into Eqn.7, we get the 
Lagrangian in the balanced case where the symmetry is enlarged to U(1)l x 17(1) g x Zi- 

III. GROUND STATE WAVEFUNCTION 

In this section, we derive the ground state wavefunction from the formalism reviewed in the last section. From 
Eqn.7, we can find the corresponding Hamiltonian in the charge sector: 

n c = \ £[4n + (-$)(v + (g) + ^ JL)n+($) + ^M-flMfl] (9) 

q 

where U+(q) = Sp+/2 and [0+(q), n+(^)] = iK5(q + (f). 

Representing 6+(q) by — ; m the long wavelength limit, neglecting the Coulomb interaction V+(q) <~ 1/q 

which is less singular than we find the ground state wavefunction in the charge sector: 

i o_ 

* b c0 [Il + (q)} = cxp [- J2-S P+ (-q)Sp + (q)} (10) 
1 1 q 

In the following, we use x to stand for the complex coordinate x + iy. Using 5p + (x) =^8{x — Zi) + ^5(x — Wi) — p 
and transforming to the position space, it is simply the modulus of the (111) wavefunction in the balanced case 
*c0 = l*ni I- This is the wavefunction in the bosonic picture. In order to get the wavefunction in the original 
fermionic picture, we need to perform the inverse of the SGT in Eqn.4 on the bosonic wavefunction. In the first 
quantization form, the inverse is: 

1L \z l - Zj\ LL \w t - wj\ ~ AJ - \Zi - Wj\ 

Performing the inverse transformation on the modulus leads to the (111) wavefunction in the fermionic coordinates: 

*co = E>b*d) = *iii («,«>) (12) 

In contrast to the SLQH, there is also an additional pseudo-spin sector in the BLQH which contains the most 
interesting physics. From Eqn.7, we can find the corresponding Hamiltonian in the spin sector is: 
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n s = ^5>f- n -(-fl n -(fl + ?*g-6_(-w_®] (is) 

where n_(g) = 6p~/2 and [0_(g),II_(<f)] = ihS(q + if); ps = p/m is the spin stiffness [21]. At small q, V-(q) = 
a — bq + cq 2 [21] where a ~ eP, 6 ~ d 2 and c remains a constant at small distances [12,16]. It is important to stress 
that this form of V-(q) has the shape displayed in the Fig.l, it not only has a phonon part near q = 0, also has a 
roton part near q = qo ~ 1/Zs where Ib is the magnetic length. 

Representing 9-(q) by — i dn , we find the ground state wavefunction in the spin sector: 

tt s0 [n_( g -)] =exp[-l]T ^- (g)/p % _(-g)gp_( g -)] (14) 

2 9 9 

It is easy to see that the above equation make senses only when V— (q) is positive for all q. 
At d = 0, a = b = 0, V_(q) = cq 2 , then Eqn.14 becomes: 

* s0 [n_(<7);cZ = 0] = exp[- t;^2V c/ p E Sp-(-q)5p- (q)\ = const. (15) 

9 

At any finite distance, as long as d < d c \ ~ Ib in the Fig.l, so the roton has a large gap, we can neglect the 
contributions from the roton part in Fig.l and only focus on the phonon contributions. In the long wavelength limit 
q <C qo ~ 1/^s, U_(q) — > a ~ d 2 . Using Sp-(x) = J2$( x ~ z i) ~~ J2$( x ~ w i) an d transforming to the coordinate 
space, it is 

V s0 (z, w) = exp[-2 v /T7^( V d - V - ,+V, ,)] (16) 

^— ' — zA ^\Zi — wA ^ \Wi — Wj\ 

Obviously, the above equation only holds at small distance d < d c \ and in the long distance limit \zi — Wj\ 3> Ib 
[22]. 

The total wavefunction is: 

Vo(z,w) = ^iu(z,w)^ s0 (z,w) (17) 

It is easy to see that the total wavefunction coincides with the (111) wavefunction only in d — > limit. At any 
finite d, it has an extra factor from the gapless spin sector ^ s0 (z 7 w). Note that this extra spin factor Eqn.16 is not 
in the LLL, this should not be too worrisome, because similar to the meron wavefunction listed in Eqn.l, modulus of 
the coordinates could appear in the long distance limit where Eqn.16 hold. 



IV. QUASI-HOLE WAVE FUNCTIONS 

By inserting one static vortex at the origin in layer 1 by setting J^ v = Jq V = S(x) or layer 2 by setting Jq~" = 
—Jq V = S(x), from the dual action Eqn.8, we will first try to derive the quasi-hole wavefunction and compare this 
quasi-hole wavefunction with the known "quasi-hole " wavefunction and "meron" wavefunction written down in [13]. 

In order to derive the quasi-hole wavefunction, we have to resort to the dual action Eqn.8 in the balanced case 
where the last term vanishes. Setting the two sources Af^ = 0, in the Coulomb gauge V • 6„ = 0, Eqn.8 becomes: 

777 777 i — * — * 

L d = -i2irb+e a0 d a b+ + i7rb+jv + + inb+J? + ^M) 2 + ^=(^) 2 + ^(V x b+)V+(q)(V x 6+) 

777 777 ' — * — * 

+ mb^r ~ +inb~r a - + Y-id^) 2 + Y -(dob-) 2 + -(V x b-)V.(q)(V x b~) (18) 

Note the absence of CS term in the spin sector. 

We only consider static vortices, so = J^T = 0. Integrating out 6g and 6g and transforming into the coordinate 
space lead to: 

Cd = ~h£n~ i7rJ ° + ~ 2 ^, d a b+)ln\x - y\{-KJ v Q + - 2ire a , p d a b+) + ^=(d b+) 2 + \{V x b+)V+(q)(V x b+) 
+ ^(5o&«) 2 - ^(n.r -)ln\x y\(nr -) + ±(V x &-)V_($)(V x 6") (19) 
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The corresponding Hamiltonian is: 



U = l^1 2 + (-q)e + (q) + \^jv+(-q) ^ + {-q))^J v + (q) 4ttII + ($)) + 2Il + (-q)V + (<?)Il + (q) 

Z lit Zi ill Q 



Z iJt Zi Til Q Z 

where [6±(q), Il±(?)} = iK5(q + 

From the above Hamiltonian in the bosonic representation, we can see the charge sector and spin sector remain 
" decoupled " [20]. Due to the absence of the C-S term in the spin sector, the inserted vortex only shifts the total 
density variable in the charge sector, but does not couple to the relative density in the spin sector, so the Hamiltonian 
in the spin sector remains the same as the ground state one Eqn.13, the corresponding wavefunction remains the same 
as the ground state one in the spin sector Eqn.14. All the effects of the inserted vortex are encoded in the charge 
sector. Again, neglecting the Coulomb interaction V + (q) in the long wavelength limit, we find that the wavefunction 
in the charge sector is: 

< qh = exp[i J2(\jo + (-<l) S P+ (-q))(-^)(^r + (q) Sp + (q))} (21) 

Transforming to the coordinate space and setting Jq + (x) = Jq~ (x) = —5(x) [20] lead to: 

1 f 1 1 

^c q h = ex P [- / dxdy(-S(x) + (^5{x - z x ) + 5{x - Wi ) - p))ln\x - y\{-6{y) + - z t ) + 5{y- w t ) - p))] 

i i 
Ni N 2 Ni N 2 Ni N 2 

= n n Ki* n n \« n \« -^u k - m) 

i i i—1 j — 1 i<j i<j 

The SGT for the quasi-hole could be different from that for the groundstate. If one inserts a vortex at the origin at 
the layer 1 in the boson Lagrangian Eqn.5, in order to recover the original electronic Hamiltonian Eqn.3, Uq in Eqn.ll 
is needed to remove the C-S term, an additional SGT U v \ is needed to remove the effects of the inserted vortex. In 
the first quantization, it is easy to show that [23] : 

Uvl =e i ^> r9 ra =e E, 1 ^ = Y[IL (23) 

The total SGT for the quasi-hole at the layer 1 is U q h = UqU v i. 

Performing the SGT on Eqn.22, we get the quasi-hole wavefunction, 

Ni 

^ qh {z, w) = (Q Zi ) J] l^l 1 / 2 ^^, w) (24) 

% i 

where ^o(z, w) is the ground state wavefunction Eqn.17 and A^i = N2 = N in the balanced case. Note that there is 
no singularity at the origin. 

Note that even in the d — ► limit, Eqn.24 differs from both the " quasi-hole " and the " meron " wavefunction 
listed in Eqn.l. This should not cause any problem. It is known that (111) wavefunction is the exact wave function 
in the d — > o limit. But both both the " quasi-hole " and the " meron " wavefunction listed in Eqn.l make sense 
only at finite distances. It is known that the lowest energy excitation at d = is a skymion carrying charge 1, while 
the " quasi-hole " and the " meron " carry total charge 1/2, so they are not valid excitations anymore at d = 0. We 
conclude the quasi-hole Eqn.24 make sense only at finite d. It is not interesting to take d — > limit to this equation 
[24]. 

Compared to the " quasi-hole " and the " meron " wavefunction listed in Eqn.l, we can see that there are two 
modifications in the Eqn.24: (1) The ground state wavefunction is the correct one Eqn.17 instead of the (111) 
wavefunction. (2) The prcfactor is different from both the "quasi-hole " and the " meron ". We expect this prefactor 
takes care of the strong interlayer correlations. All the wavefunctions in Eqn.l arc built upon the (111) wavefunction. 
As suggested in [13] and explicitly shown in this paper, (111) wavefunction is not even qualitatively correct at any 
finite d. As shown in [13], although prefactors are not important in SLQH due to the gap in the bulk, they maybe 
crucial in BLQH due to the gapless mode in the interlayer correlations. The two factors maybe responsible for the 
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" quasi-hole " 's charge distribution spreading over the whole system and its energy diverges linearly with the area 
of the system. Although, the "meron " wavefunction's energy is only logarithmically divergent, it ignores the strong 
interlayer correlations, so it is not a good trial wavefunction either. We propose that the quasi-hole wave function 
Eqn.24 not only has logarithmically divergent energy, well localized charge distribution, but also correct interlayer 
correlations. 



V. A PAIR OF QUASIHOLE EXCITATIONS WITH CHARGE 1 

Now we put two vortices into the BLQH system. One is on the top layer at zo and the other is on the bottom 
layer at wq. The only change from the quasi-hole calculation in the last section is Jq^ = 5(x — zq) ± 8(x — Wo) [20]. 
Again, due to the lack of the CS term in the spin sector, the Hamiltonian in the spin sector remains the same as the 
ground state one Eqn.13, so the wavefunction is not affected at all by the insertion of the two vortices and remains 
the same as the ground state one in the spin sector Eqn.14. All the effects of the inserted two vortices are encoded in 
the charge sector. The wavefunction in the charge sector in the bosonic picture is: 

^ b cl = exp(- / dxdy[-(S(x - z ) + S(x - w )) + S(x - Zi) + 5(x - Wi) - p)]ln\x - y\ 
[{\{${y- zo) + S(y- w )) + (J2s(y- Zl ) + 5{y- w t ) - p))} 

i 

ATi JVi N 2 N 2 

= Y[ \zi - z |2 Y[ \zi - w \* J| \wi - zo| 3 Y[ \wi - w | 3 (25) 

i i i i 

It is easy to see that in the bosonic picture, the above wavefunction in the charge sector is symmetric under Zi <-» 
or z <-» w separately. This is under expectation, because the two layers are completely symmetric in the charge 
sector. 

Just like deriving U v \ for the quasi-hole, we can get U V 2 for a pair of vortices inserted at z in top layer and wo at 
the bottom layer: 

\Zi - Z \ \Wi - Wq\ 

% 

The total SGT for the meron pair is U pa ir = UqU V 2- 

Performing the SGT on Eqn.25, we get the wavefunction for a pair of quasiholc: 

V pair (z, w;z ,w Q ) = (Wzi - z )( Wi - w ) TT | ^ ~ ~ ^ \V 2 * Q (z, w) (27) 

i (Zi - z )(w l - z ) 

where ^o(z,w) is the ground state wavefunction Eqn.17 and N± = N 2 — N in the balanced case. Note that the pair 
wavefunction is not symmetric under zi <-» Wi or z <-> w separately anymore. This is because the SGT U v2 Eqn.26 
is not. Of course, it is still symmetric under Zi <-> Wi, z <-» w simultaneously. 

If we insert the two vortices at the same point, namely, putting z = w = in the above equation, as expected, 
we get: 

^ P air(z, W, 0) = (JJ ZiWi)Vo(z, w) (28) 
i 

This corresponds to insert a single vortex through the two layers. In contrast to the quasi-hole excitation Eqn.24, 
Eqn.28 carries charge 1 and remains a valid wavefunction even at d = 0. Indeed, in the d — > limit, it recovers the 
" meron pair " wavefunction listed in Eqn.l. If one splits the whole vortex, it will evolve into the pair wavefunction 
Eqn.27. 

The "pair meron" wavefunction built on (111) is listed in 1( essentially Eqn.(llO) in Ref.10 ). It was shown that its 
energy E pair ~ \z a — w \ 2 instead of logarithmically as naively expected, because the charges are extended between 
Zo and w a - Similar to the quasi-hole wavefunction Eqn.24, there are two modifications in the Eqn.27: (1) The ground 
state wavefunction is the correct one Eqn.17 (2) The prefactor is different. We expect this prefactor takes care of the 
strong interlayer correlations between the two vortices, the pair wave function not only has logarithmically divergent 
energy, well localized charge distribution, but also correct interlayer correlations. 
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VI. INSTABILITY IN THE WAVEFUNCTION AS THE DISTANCE INCREASES 



When the distance is sufficiently small, the BLQH is in the ESF phase, we expect the ground state, quasi-hole 
and pair wavefunctions Eqns. 17,24,27 only hold in the ESF phase. When the distance becomes sufficiently large, 
the two layers become two weakly coupled Fermi liquid ( FL ) layers. All these wavefunctions completely break 
down. New set of wavefunctions are needed. Although the ESF phase and FL phase at the two extreme distances 
are well established, the picture of how the ESF phase evolves into the two weakly-coupled FL states was not clear, 
namely, the nature of the intermediate phase at d c \ < d < d C 2 was still under debate. Recently, starting from the 
well studied excitonic superfluid (ESF) state, as distance increases, one of the authors found [14] that the instability 
driven by magneto-roton minimum collapsing at a finite wavevector in the pseudo-spin channel leads to the formation 
of a pseudo-spin density wave (PSDW) at some intermediate distances. He constructed a quantum Ginsburg-Landau 
theory to study the transition from the ESF to the PSDW and analyze in detail the properties of the PSDW. He 
showed that a square lattice is the favorite lattice. 

As shown in [14] , it is the original instability in V- (q) = a — bq + cq 2 which leads to the magneto-roton minimum 
in the Fig. la. By looking at the two conditions V-(q)\ q=qo = and dV d ^ \ q = qo = 0, it is easy to sec that V-{q) 

indeed has the shape shown in Fig. lb. When b ~ d 2 < b c = ~ d, the minimum of V-(q) at q = qo = \fafc ~ d 

has a gap, the system is in the ESF state, this is always the case when the distance d is sufficiently small. However, 
when b — b c , the minimum collapses and S(q) diverges at q = q , which signifies the instability of the ESF to the 
pseudo-spin density wave (PSDW) formation. When b ~ d 2 > b c — 2 v / ac <~ d, the minimum drops to negative, the 
system gets to the PSDW state, this is always the case when the distance d is sufficiently large. 




ESF dcl PSDW dc2 FL 

<\|/>^0 <\|/> = 

<n7r>=0 <n 7 ->^0 

<j G 



Fig.l: The zero temperature phase diagram in the balanced case as the distance between the two layers increases. ESF where 
< 4> 0, < riQ >= stands for excitonic superfluid, PSDW where < tp >= 0, < >^ stands for pseudo-spin density wave phase, FL 
stands for Fermi Liquid, (a) Energy dispersion relation ui(q) in these phases, (b) Ve(<?) in these phases. The cross in the PSDW means 
the negative minimum value of Ve(q) is replaced by the PSDW. The two order parameters were defined in [16]. In reality, the instability 
happens before the minimum collapses. 

We can easily see the instability from the ground state wavefunction Eqn.17. The only distance dependence in the 
charge sector is encoded in the V + {q) in the bosonic Hamiltonian in the charge sector Eqn.9, but this dependence is 
ignored in the (111) wavefunction in the charge sector. The d dependence in V+ is smooth anyway. Essentially all 
the distance dependence is encoded in the ground state wavefunction in the spin sector Eqn.14. As can be seen from 
Fig.l, the instability happens at q = qo where V— becomes negative, but the spin stiffness pe remains non- critical 
through the ESF to PSDW transition. As d — > d~ 1: the sum over q in Eqn.14 becomes dominated by the regime 
q <~ qo- When d c \ < d < d C 2, Eqn.14 breaks down. A new wavefunction to describe the translational symmetry 
breaking PSDW state is needed. Some trial wavefunctions are proposed in [25]. It would be interesting to derive the 
new wavefunction of the PSDW state from the CB theory. 

VII. CONCLUSIONS 

BLQH differs from the SLQH by its symmetry breaking ground state and associated neutral gapless mode in the 
pseudo-spin sector. Due to the gapless mode in the bulk, the groundstate wavefunctions could be considerably different 
from the well known (111) wave function [12,13]. The low energy excited states could also be sensitive to details such 
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as normalization factors. One important problem is to find good trial wavefunctions for the ground state and low 
energy excited states. We investigated this important open problem from the CB theory developed previously in 
[14,16] to study BLQH systematically. We derived the ground state, quasi-hole and a pair of quasi-hole wavefunctions 
from CB theory and its dual action by the following procedures: We first performed the singular gauge transformation 
Eqn.4 to transform a fermionic problem into a bosonic problem, then found that all the wavefunctions in the bosonic 
picture are always the product of two parts, one part in the charge sector and the other in the spin sector. All the 
distance dependence are encoded in the spin part, while all the excitations only happen in the charge sector. After 
transforming back to the original electron picture by proper inverse SGT's, we get the final wavefunctions in the 
electron coordinates. We found that the inverse SGT's are different in ground state, meron and a pair of merons. By 
considering the differences carefully, we derived all these wavefunctions in the original electronic picture a systematical 
way. 

At any finite d, the ground state wavefunction in the charge sector is the same as the (111) wavefunction, while that 
in the spin sector is highly nontrivial due to the gapless mode. So the total groundstate wavefunction differs from the 
well known (111) wavefunction at any finite d. In the bosonic picture, when inserting vortices in the ground state, 
the spin part remains the same due to the lack of CS term in this sector, while the charge part changes accordingly. 
However, due to the insertion of vortices, in order to recover the original electronic problem, the inverse SGT differs 
from that in the ground state. After transforming back to the original electron problem by the inverse SGT's, we 
showed that the quasi-hole and a pair of quasi-holes wavefunctions contain non-trivial normalization factors as shown 
in 24,27. We expect that the quasi-hole and pair wave function not only has logarithmically divergent energy, well 
localized charge distribution, but also correct interlayer correlations. It is important to test these trial wavefunctions 
by QMC simulations performed in [13] for the states listed in Eqn. 1. We also investigated the instability encoded 
in the spin sector which leads to the PSDW solid formation proposed in [14,16]. Because the CB field theory has 
been used to describe the tri-layer quantum Hall systems very successfully, the analysis in this paper can be easily 
extended to derive the wavefunctions in the TLQH [17]. 

It is well known that CB approach is not a Lowest Landau Level ( LLL ) approach [12,16], it is very difficult to 
incorporate the LLL projection into the CB approach. This may be partially responsible for the spin part of the 
ground state wavefunction Eqn. 16 not in the LLL level. But as explained below Eqn. 17, this should not be too 
worrisome, because Eqn. 16 works only in long distance anyway. As shown in [13] and listed in Eqn.l, in the long 
wavelength limit, the meron wavefunction's normalization factor contains modulus which is not in the LLL cither. 
How to getting precise short distance behaviors of these wavefunctions from the CB theory remains an open problem. 

J. Ye thanks S. C. Zhang for helpful discussions. He also thanks Yu Lu and Xiang Tao for the hospitality during 
his visit to the Institute for Theoretical Physics of China. 
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